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ABSTRACT 



We examine a sample of 223 F, G and early K metal-poor subdwarfs ([m/H] < 
— 1) with high proper motions (/i > 0.2"/ year) at the distances of up to 250 
pc from the Sun. By means of our own speckle interferometric observations 
conducted on the 6 m BTA telescope of the Special Astrophysical Observatory 
of the Russian Academy of Sciences, and the spectroscopic and visual data taken 
from the literature, we determine the frequency of binary and multiple systems in 
this sample. The ratio of single, binary, triple and quadruple systems among 221 
primary components of the sample is 147:64:9:1. We show that the distribution of 
orbital periods of binary and multiple subdwarfs is asymmetric in the range of up 
to P = 10 10 days, and has a maximum at P = 10 2 — 10 3 days, what differs from 



199l|). 



the distribution, obtained for the thin disc G dwarfs ( iDuquennoy fc Mayor! !" 
We estimated the number of undetected companions in our sample. Comparing 
the frequency of binary subdwarfs in the field and in the globular clusters, we 
show that the process of halo field star formation by the means of destruction of 
globular clusters is very unlikely in our Galaxy. We discuss the multiplicity of 
old metal-poor stars in nearby stellar streams. 



Subject headings: stars: binaries including multiple: close — stars: binaries: 
spectroscopic — stars: binaries: visual — stars: Population II — stars: subdwarfs — 
Galaxy: halo 
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1. Introduction 

The study of metal-poor stars ([Fe/H] < —1) makes it possible to shed light on 
numerous problems of modern astrophysics, such as heavy elements production in supernova 
explosions, the metallicity distribution function of the stellar halo, the initial mass function, 
the nature of the Big Bang and the first generation, or Population III, stars, etc. An 
important place among these fundamental problems is occupied by the questions of the 
origin, chemical and dynamical evolution of our Galaxy. The oldest stars with the masses 
of M ^ 0.8 M & are unevolved. Therefore, the abundance of chemical elements in their 
atmospheres reproduces the composition of prestellar matter. Additional information on 
spatial motions of these stars preserves the possibility to reconstruct the way the Milky 
Way formed. 

The orbital elements of binary and multiple stellar systems are an important tool 
for studying prestellar matter. In single low mass stars, the mass is the only parameter 
conserved since the time of star formation. Binary and multiple systems bear three more 
conserved va lues: the ang ular momentum, the eccentricity and the mass ratio of their 



components (ILarson 



2001 



J) in case of detached systems. Therefore, binary and multiple 
stars carry more information on the process of star formation than single stars. The 
study of binary and multiple metal-poor systems enables us to impose certain restrictions 
on the physical conditions in prestellar matter at the time of the genesis of our Galaxy. 
Metal-poor stars are common in the globular clusters, galactic halo and i n the galactic 



Eggen 



1996a 



field, where an existence of the so-called stellar streams was revealed (e.g. 
The multiplicity and the orbital parameters of binary and multiple stars in these streams 
may also provide additional information on the nature of the stream's progenitor and its 
dynamical evolution. 

The problem of stellar multiplicity was widely discussed in the literature, however, 
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it mostly concerned the thin disc stars with so 



1991 



Fischer fc Marcy 



1992 



Halbwachs et al. 



ar-li 



ce metallicities (IDuquennoy fc Mayor 



20031 ). Metal-poor stars were studied 



much less, as the i r occ urrence in the solar neighbourhood is less than 1%, according to 



Nordstrom et al. 



( 120041 ) catalog. Early works, addressing the rate of binary systems among 



Population II stars, sh owed that this va^ 



for P opulation I stars (lAbt fc Levi 



1987). In subsequent works (IPreston fc Sneden 



1969 



ue is small, as compared to the analogous value 



Crampt on fc Hartwick 



2000 



1972 



Goldberg et al 



Abt fe Willmarth 



2002; 



20021 ) it was concluded that these values are indistinguishable (idem 



Abt 



Latham et al. 



2008!). A 



l ong-term spectrosc opic monitoring of a 



DOUt 



( jCarney et al. 



1994! (hereinafter CLLA). 



2001 



500 nearby stars with high proper moti ons 



Goldberg et al. 



2002 



Latham et al. 



20021 ) has 



played an important role in the study of the multiplicity of metal-poor stars. Spectroscopic 
studies cover the systems with relatively short orbital periods (P < 10 years). 

The study of long-period co uples with common proper motion components 



(jZapatero Osorio fc Martin 



20041 ) confirms the hypothesis of an equal frequency of binary 



Allen. Poveda fc Herrera 



2000). 



systems among the old and young stellar populations (idem 
Meanwhile, an 'intermediate' period range of P ~ 10 — 1 000 years, which corresponds to 
the semi-major orbital axes of a ~ 10 — 100 AU in the solar neighbourhood, remains poorly 
understood to date. This range can be studied with the use of adaptive optics, speckle 
interferometry and long baseline interferometry. The scarce Population II stars observations, 
made by means of the interferometric techniques, were ra n either for the brightest star s 



(|Lu et al. 



19871 ) or with relatively low angular resolution ( jZinnecker. Kohler fc Jahreifi 



20041 ) . Notwithstanding the empirical data available to date, the number of known binary 



and multiple systems with metal-poor components remains small. 

In order to enlarge the database of binary and multiple Population II stars, to define 
their orbital parameters and the properties of their components, we conducted speckle 
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interferometric observations of 223 metal-poor subdwarfs with high proper motions loc ated 



in th e solar neighbourhood ( IRastegaev. Balega &: Malogolovets 



2007 



Rastegaev et al. 



20081 ). The observations were made with the diffraction-limited resolution of the 6 m Big 



Telescope Alt-azimuthal (BTA) of the Special Astrophysical Observatory of the Russian 
Academy of Sciences (0.023" at the wavelength of 550 nm). In the present work, we analyse 
the multiplicity and orbital periods distribution for binary and multiple stars. We made our 
analysis based on own observations and the data adopted from other authors. Additionally, 
an attempt was made to examine the ratio of binary and multiple stars in the streams of 
old metal-poor stars located in the solar neighbourhood. 



2. Sample 



For the observations with high angular resoluti on, we compiled a sam ple of 223 field 



subdwarfs of the F, G and early K spectral classes (IRastegaev et al 



20071 ) from the CLLA 



catalog. The CLLA presents a spectroscopically stu died sample of the A-K spectral 



types dwarfs from the Lowell Proper Motion Survey (IGiclas. Burnham fc Thomas 



1971 



19781 ). which mainly includes the stars from the Northern Hemisphere with proper motions 



exceeding 0.26" per annum and brighter than 16 m in the B band. 



We selected 223 stars from the CLLA using the following criteria: 



metallicity [m/H] < — 1, 



declination 5 > —10°, 



apparent magnitude my < 12 r 



The last criteri on was determined b y the limiting stellar magnitude of our speckle 



interferometer (IMaximov et al. 



2003), which was about 13 m . No restrictions were applied 
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on the heliocentric distances of these stars, evenly distributed on the celestial sphere. The 
maximum distance to the sample objects is 250 pc. The median heliocentric distance of the 
selected stars is approximately 100 pc. This allows us to take advantage of high angular 
resolution in order to detect new systems, since at such distances the pairs with semi-major 
orbital axes from 10 to 100 AU are hard to detect both spectroscopically and visually. 

Using the two following criteria, fi > 0.26 "/year and [m/H] < —1, we tried to rule out 
the thick and thin disc stars. In case of the thick disc though, these limi ts are not string ent 



and s ome objects may belong to the metal-weak tail of the thick disc ( lArifyanto et al 



2005|). 



To separate the halo sta rs from the disc stars in our s ample, we used a formal method 



described in the appendix of 



Grether fc Lineweaverl (120071 ) . The equations establishing the 



probability that a star belongs to the thin disc (fthin)j the thick disc (Pthick) or the halo 
(-Phaio) are 

P 1 P 2 P 3 

-Pthin — fl— , -Pthick = f2~p > -Phalo = f'S~j^i 



where 



Ci exp 



u 2 (v-(v)y 



w 2 



2< 



24 



2 < 



([Fe/H] - ([Fe/RW 

2ct 2 

[F e /H]j 



1 (thin disc), 2 (thick disc), 3 (halo). 



i r 



The input parameters for these equations are adopted from iRobin et al.l (120031 ) . and 
presented in Table [TJ We rejected three out of 221 systems: two double stars G99-48 and 
G166-45, and one single BD — 1°1792, because for these systems there are no heliocentric 
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distances or UVW components in the CLLA. We referred the objects with Phaio > 0.5 
to the halo stars. There is a total of 148 of such stars in our sample. The remaining 70 
objects belong to the thick disc. None of the stars from the sample have Pthin exceeding 
0.01. Average metallicity and space velocity vector components for halo stars in our sample 
are: {{[Fe/R] halo ), (U halo ), (V halo ), (W halo )) = (-1.9 ± 0.5, -16 ± 152 km/s, -180 ± 81 km/s, 
—1 ± 78 km/s), where the errors are standard deviations. In a similar way, these values for 
the thick disc stars are: (([Fe/H] tWcfe ), {U thick ), (V thick ), (W t hick)) = (-1-2 ± 0.2, -12 ± 83 
km/s, —92 ± 59 km/s, 3 ± 50 km/s). However, the system of equations ([1]) might be unable 
to accurately describe the situation in the areas where the star parameters of different 
Populations overlap. 

Fig. [1] shows the distribution of the sample stars on the graph metallicity versus 
V^-component of spatial velocity, data taken from the CLLA. The halo stars with a 
substantial dispersion of spatial velocities and with low metallicities are located in the 
central and left parts of the figure. At the top rig ht of the figure, along with the halo stars, 



there are stars of the thick disc's metal-weak tail (lArifyanto et al. 



2005|). 



The sample consists of main sequ ence stars and seven blue stragglers (jCarney et al. 



2001 



Carney. Latham fc Laird 



2005b! ) . — a continuation of the main sequence in the area 



of hotter and bluer stars, as compared to the turnoff stars. Fig. [2] represents the location of 



the sample stars in the Hertzsprung-Russel 



M13 globular cluster with [Fe/H] = —1.61 (IGrundahl. VandenBerg fe Andersen 



diagram compared to the old popu 



at ion of the 



19981 ) 



Fig. [3] represents a comparison between the trigonometric and photometric distances 
for the sample stars. The parallax es for 135 sample stars were taken from the HIPPARCOS 



satellite data ( van Leeuwen 



20071 ). The corresponding photometric distances for 133 stars 
are retrieved from the CLLA. We can see that the distances, obtained with HIPPARCOS, 
are systematically larger than the corresponding distances from the CLLA. This difference 
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can be explained by some unaccounted components, the presence of which may lead to 
an underestimation of photometric distances. Another reason for this discrepancy is the 
photometric distance calibration adopted in the CLLA. 



Most o 



( Allen et al 



the s t ars in our sample were examined for common proper motion components 



2000; 



Zapatero Osorio fc Martin 



20041 ) . 



3. Observations and results 



The speckle interferometric obse rvations of 223 sample stars were carried out in 



2006-2007 on the 6 m BTA telescope (IRastegaev et al. 



2007 



20081 ) which diffraction-limited 



resolution is 0.023" for A = 550 nm and 0.033" for A = 800 nm. Most of the observations 



were carried out using the system (IMaksimov et al. 



20091 ) based on a 512x512 EMCCD (a 



CCD featuring on-chip multiplication gain) with high quantum efficiency and linearity. This 
system allowed us to detect objects with magnitude differences between the components 
of up to Am = 5 m . Taking into account the limiting stellar magnitude of our sample 
(my < 12 m ), detected secondary component can be as faint as 17 m . The 4.4" field of view 
of our system allows detection of secondary components at a separation of 3" from the 
primary star. The speckle interferograms were recorded using five filters: 545/30, 550/20, 
600/40, 800/100 and 800/110 nm (the first number indicates the central wavelength of the 
filter, the second — the half-width of the filter's bandwidth) with the exposures of 5 to 20 
ms. For each object, we accumulated from 500 to 2 000 short exposure images depending 
on weather conditions. The observations were made with an average seeing o f 1.5". The 



20021 ) may be as 



accuracy of our speckle interferogram processing method (iBalega et al. 
good as 0.02 m , 0.001", and 0.1° for the component magnitude difference, angular separation 
and position angle respectively. 
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For 19 stars in our sample we observed the speckle interferometric companions. Sixteen 
companions were resolved astrometrically for the first time. We discovered 5 new binary 
systems (G191-55, G114-25, G142-44, G28-43, G130-7), 3 triple systems (G87-47, Glll-38, 
G190-10) and one quadruple system, G89-14. The position parameters and magnitude 
differences between the speckle interferometric components are listed in tables 1 and 2 of 



Rastegaev et al. 



(120081 ) . 



4. Sample completeness 

To be able to determine the completeness of our sample, we used the star count 
method. From all known subdwarfs within 25 pc from the Sun, we selecte d 5 objects with 



abso 



Fuchs fc Jahreifi 



ute magnitudes My < 8 m and metallicities [Fe/H] < — 1 (see tab. 1 in 
19981 ). The star GJ 1064 A, with the metallicity of —1 dex, was as well included in our 
sample. An extrapolation to the volume of our sample (with the radius of ~ 250 pc) 
increases the number of such objects to ~ 5 000. Obviously, when we consider the frequency 
of stars in different volumes, we cannot depart from the uniform distribution of stars in 
space. We have to bear in mind the structure of the Galaxy and the data available on the 
star distribution in its various subsystems. However, as we show in the Appendix, for the 
stars in our sample located within 250 pc from the Sun, the structure of our Galaxy can 
be neglected. Therefore, the number of objects in our sample constitutes about 5% from 
the total number of stars in question located in the examined region of space. We have to 
mention that while examining the stars in different volumes, we were taking account of the 
primary components only in case of spectroscopic and speckle interferometric pairs, and of 
both components in the case of visual and common proper motion pairs. 
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5. Undetected companions 

Capability to detect a binary system is determined by both the observations method 
used and by the physical characteristics of the system itself. As the sample objects were 
examined using three different methods, — spectroscopic, speckle interferometric and visual, 
we have to analyse the number of systems unaccounted for by each of these methods. To do 
that, we conditionally divided the spectroscopic and astrometric pairs by the values of their 
semi-major axis at a < 10 AU and a > 10 AU, and their orbital periods at P < 10 000 days 
and P > 10 000 days, respectively. 

For spectroscopic pairs, we hypothesized that the component mass ratio is uniformly 
distributed from to 1, then took the deduced period distribution for the studied stars in 
the r ange of to 10 000 days (see belo w), and made our calculations using the formulae 



from 



Mazeh. Latham fc Stefaniki (119961 ). Our conclusion is that the probability of not 



detecting a binary system approximately equals to 20%. 

As for the visual and interferometric binaries, the secondary can be detected if 
the following two conditions are satisfied. Firstly, the magnitude difference between 
the components should not exceed a certain critical value, determined by the detector's 
dynamical range and the spectral range used. Secondly, in case of speckle interferometry, 
the angular separation between the components should be larger than the telescope's 
diffraction limit and smaller than the detector's field of view, or, in case of visual studies, 
the angular separation should be smaller than the certain restrictions imposed. 

In order to estimate the number of unresolved systems due to the effect of ellipse 
projection of the true orbit on the picture plane, and due to the secondary component's 
orbital phase at the moment of observation, we used the Monte Carlo simulation. We 
modelled a sample of stellar orbits, evenly distributed in space, with heliocentric distances 
from 25 to 250 pc, uniformly distributed eccentricities (from to 0.9) and arguments of 
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periapsis (from to 360 °). Inclinations of the true orbits to the picture plane obey the 
sin i law. 

The result ing orbits were projected onto the picture plane using the following formulae 



(ICouteau 



19811 ) 



a (1 — e 2 ) cos (v + oj) 
^ 1 + e cos v cos (9 — Q) ' 

tg (9 - n) = tg (v + u) cos i, (2) 
v ll + e u 
tg 2 = VT^ tg 2' 
M = u — e sin u, 



where p is the angular separation between the components in the picture plane, a is the 
orbital semi- major axis (in arcseconds), e is eccentricity, v is the true anomaly, u is the 
eccentric anomaly, M is the mean anomaly, u is the argument of periapsis, 9 is the position 
angle of a secondary component, Q is the longitude of an ascending node and i is the 
inclination angle between orbital and picture planes. 

Further, a random orbital position of the satellite was predetermined for each simulated 
system. At this point we accounted for the fact that a companion spends more time 
at apoastron than at periastron. Then the projection distances p from the satellite to 
the primary star were counted. A relative number of projection distances greater than 
0.033"(the diffraction limit of the 6 m telescope in a 800/100 filter) is, in fact, the probability 
value of a system detection using the speckle interferometric and visual methods. The 
results of modelling show that for the systems with an orbital semi-major axis exceeding 
a = 10 AU, such probability is close to 1 and an influence of the components' geometry can 
be neglected. This probability weakly depends on the choice of a plausible distribution of 
semi-major axis for a > 10 AU. 

Let's consider now the incompleteness of detection caused by the magnitude difference 
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between the components. An average mass of a primary star in our sample is 0.65 M & 
and a standard deviation is a M ~ 0.05 M Q . The maximum magnitude difference between 
the components, detec t ed by speckle interferometric an d visual methods, is about 5 m 



( iMaksimov et al. 



2009 



Zapatero Osorio fc Martin 



2004 ). For speckl e interferomet r y, we 



Baraffe et al 



(1997) 



used a more conservative estimate of 4 m . According to the model of 
for [Fe/H] = —1, the minimal mass of a secondary, which is 4 mag fainter than a primary 
of an average mass, equals M m i n = 0.20 M & in the I band. This band roughly corresponds 
to our 800/100 and 800/110 filters. To compute the number of undetected companions, we 
have to know the function of mass ratio distribution. We used a uniform distribution. If 
the mass ratio for wide (speckle interferometric, visual and CPM) pairs obeys a uniform q 
distribution, then from the following ratio 



M 2 =M min 

J dM 2 

M 2 =0.08 
M 2 =0.8 

/ dM 2 

M 2 =M mi „ 



(3) 



we find the quantity of secondary components with masses in the range from 0.08 M & (mass 
of a brown dwarf) to M min per one secondary component with the mass in the range from 
M m i n to 0.8 M (a maximal mass of the stars in our sample). This quantity equals 0.2. 
Therefore, a fifth of discovered speckle interferometric companions remains undetected at 
the given distribution of f(q) with the use of 800/100 (or 800/110) filter. For the 550/20 
and 545/30 filters, M min = 0.28 M and the expression ([3]) is approximately 0.38, which 
makes these filters less suitable for our task in the sense of magnitude difference (yet, more 
suitable from the viewpoint of the angular resolution). The se arch for common proper 



motio n pairs for the sample stars was conducted in the I band (IZapatero Osorio fc Martin 



20041 ) . Forty seven objects from our sample were observed speckle interferometrically, solely 
using the 550/20 or 545/30 filters, 9 objects were observed only in the 600/40 filter and one 
object (G183-9) was observed in the 550/20 and 600/40 filters. The 166 remaining stars 
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were at least once observed in the 800/100 or 800/110 filters. With this in mind, we find 
that the expression ([3]) is approximately equal to 0.25 for all filters. Therefore, according 
to these rough estimates, for 33 astrometric components in our sample (see Fig. H]) we 
have about 8 unaccounted companions or 24%; that is comparable to 20% of unaccounted 
spectral companions. Thus, we can assume that the number of unaccounted components in 
our sample does not depend on the orbital period and is 20 — 25%. 



Analogously, the IMF-like distribution f(q) ~ for 0.08 < M/M© < 0.5 ( IKroupa 



2001 ) gives us 37 unaccounted companions. In case f(q) grows towards bigger q values 



( ISoderhjelm 



20071 ). then the number of undetected astrometric components in our sample 
does not exceed 5. For our further estimates we chose the values corresponding to the 
uniform q distribution. 



6. Multiplicity of the sample 



6.1. Raw estimates 



In order to calculate the ratio of the systems of different multiplicity among the 
studied stars, we complemented the results of our speckle interferometric measurements 
by the data from spectroscopic and visual studies found in literature. The information on 
the spectroscopic companions was adopted from the publications dedicated to long-term 



spectroscopic 



Carney et al. 



monitoring of Population 



2001 



Goldberg et al. 



wide visual pairs was taken from 



the WDS catalog ( jMason et al. 



2002 



stars by Carney, Latham, Laird et al. (CLLA 



Latham et a 



Allen et al. 



feOOOj ) and 



2002 



Latha 



m 



2008). The data on 



Zapatero Osorio fc Martin 



( 12004] ) via 



20011 ) . As a result, the ratio of sing 



(S : B : T : Q) systems for the stars in our sample is 147:64:9:1 ( iRastegaev et al. 



e : binary : triple : q uadrup le 



2007 



20081 ). Therefore, at least 159 stars from 306 stars in our sample (221 main components 
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and 85 satellites) belong to binary and multiple systems. The multiplicity of our sample 



33%^ly, where multiplicity is understood as the ratio of binary and multiple systems to the 



-6%> 

total number of systems 



B + T + Q 

J systems 



S+B+T+Q 

A similar, unadjusted for unresolved c ompanions value for the thi n disc stars of spectral 



classes from F7 to G9, equals to 51:40:7:2 (IDuquennoy fc Mayor 



19911 ). with the multiplicity 



of about 50%. It is necessary to pay attention to the differences between the two compared 
samples. When our sample was formed from the stars with a certain apparent magnitude 
limit and high proper motions, the sample from Duquennoy and Mayor is only limited by 
the heliocentric distances (all their stars are located within 22 pc from the Sun). 

In terms of multiplicity, the Population II stars differ from the Population I. At least 
a third of metal-poor systems and at least half of the systems with solar like abundances 
are binary and multiple. This discrepancy can be explained by both the complexity of 
detection of low mass metal-poor spectroscopic satellites, by the selection effects, and by 
the dynamical evolution of binary and multiple stars. 



6.2. Corrected estimates 

To estimate the true multiplicity of the sample stars, we have to account for various 
selection effects that are unavoidable in astronomical observations. The five underlaying 
criteria for the choice of stars in our sample are: proper motion (fi > 0.26 "/year), 
metallicity ([m/H] < —1), magnitude (m v < 12 m ), spectral classes (F, G and early K), and 
position in the sky (all our stars are located in the Northern hemisphere). We also have to 



Everywhere further, we use a 95% confidence interval as an error for the multiplicity 
obtained by us. For this purpose we use the properties of the binomial distribution. 
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take into account the undetected companions (Section [5]). 



For the sake of simplicity we assume that the multiplicity of a stellar Population does 
not depend on proper motion and position in the sky. This gives us grounds to disregard 
the selection effects related to the proper motion of the stars and their coordinates. 

To account for the bias incurred by magn itude, we have to consider the Opik effect 



f lOpik 



1923 



Goldberg. Mazeh fc Latham 



20031 ). which applies to binaries detected within a 



magnitude-limited sample of stars. This effect is an expansion of the Malmquist effect for 
binary stars. Binaries are on the average brighter than singles. Thus, in a magnitude-limited 
sample the binaries are observed from a bigger volume in space than the single stars. To 
avoid this effect we rejected the binaries which are brighter than 12 m due to the contribution 
of the secondary component to the total luminosity. For SB1 binaries we adopted that 
the lumi nosity of the secondar y component is not less than 2 magnitudes fainter than the 



primary ( [Goldberg et al. 



20021 ). We discarded the SB1 pairs whose primary component 
was fainter than 12 m in the assumption that the secondary component is 2 magnitudes 
fainter than the primary. For SB2 pairs the luminosity of the primary component was 
determined from the t otal luminosity of the system an d the mass ra ti o of t he components 



( jGoldberg et al. 



20021 ) using evolutionary tracks from iBaraffe et al.l ( 119971 ). We thus 



rejected four binaries, one SB1 (G242-14) and 3 SB2 (G86-40, G99-48, G183-9) as the 
measured luminosities of the primaries are fainter than 12 magnitude in the l/-band. In the 
worst case for SB2 systems we have to discard 4 pairs in the assumption that for a given 
integrated magnitude the luminosities of primary and secondary components are equal. 
None of the speckle interferometric pairs have been excluded from consideration due to the 
Opik effect. The magnitude differences between the components measured by us for a given 
integrated magnitude suggest that the primary components of the speckle pairs are brighter 
than 12 m . Not a single system of higher multiplicity was dropped from the analysis as their 
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primaries are not compliant with the condition my > 12 m . 

Taking into account an adjustment for unresolved components and Opik effect, the 
multiplicity of F, G and early K subdwarfs in the solar neighbourhoo d is, according to our 



calcu lations, at least 40% and at least 60% for the thin disc G dwarfs (IDuquennoy fc Mayor 



199ll ). In both cases correction does not exceed 10%. 



6.3. Multiplicity as a function of metallicity 

To investigate the effect of metallicity on the multiplicity of stars in our sample, we 
divided it into four metallicity bins: (-3,-2.5], (-2.5,-2.0], (-2.0,-1.5], (-1.5,-1.0). 
One single star, G64-12 with [m/H] = —3.52 was rejected from the analysis. In every bin 
we evaluated the ratio of single:binary:triple:quadruple systems and multiplicity f sys tems- 
The acquired results are presented in Table [2j whereof it is clear that in the range [m/H] 
from —2.5 to —1.0 the multiplicity does not depend on metallicity and constitutes one 
third. In the range (—3, —2.5] the multiplicity is somewhat lower, but the number of 
systems in this range is smaller than in the rest of ranges, which affects the accuracy 
of determination of f systems (see Fig. |5]). It can be seen from Fig. |5] that in the first 
approximat ion the rate o f binary and multiple systems in our sample feebly depends on 



metallicity. 



Carney et al. 



(l2005al ) obtained analogous results based on a bigger sample of 
stars, studied spectroscopically. Note that the bulk of our triple stars belong to the range 
[m/H] = (-1.5,-1.0). 



6.4. Halo versus thick disc 

Using the set of equations (CQ) we separated the stars of the thick disc and the 
halo stars in our sample. The ratio single:binary:triple:quadruple systems for thick disc 
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stars is 46:19:5:0, and the multiplicity /systems = 34%^ii% • The same for the halo stars 
S:B:T:Q =101:45:4:1, and f systems = 33%+^. 

An independence of the ratio of binary and multiple systems from the conditional 
division 'thick disc-halo', as well as the consistency of f sys tems ~ 33% under metallicity 
changes (Table [2]) might testify that most of the stars in our sample belong to one and the 
same galactic subsystem, the halo. Note again that the set of equations ([I]) might not be 
operable in cases when the stellar parameter ranges of different galactic subsystems overlap. 
It is not impossible though that the Population II stars, both the halo and thick disc stars 
have similar ratios of the systems of different multiplicity. 



6.5. Multiplicity as a function of kinematics 

Table [3] represents the basic facts on the dependence of multiplicity of the sample stars 
on kinematics. We excluded three objects from consideration: two double stars G99-48 and 
G166-45, and one single BD — 1°1792, because for them there are no distances or UVW 
components in the CLLA. For single and binary and mltiple stars the average values and 
standard deviations of the space velocity vector components are: (U s ) = — 13 ± 137 km/s, 
(U b+m ) = -18± 127 km/s, (V s ) = -160±86 km/s, (V b+m ) = -135±81 km/s, (W s ) = 2±74 
km/s, (Wb + m) = — 3 ± 61 km/s. For norm of velocity vector v = y/U 2 + V 2 + W 2 : 
(v s ) = 222 ± 89 km/s, (vb+ m ) = 196 ± 80 km/s. In the first approximation we may consider 
that for our stars the multiplicity is constant and does not depend on the spatial velocity 
vector components. 

The distribution of the norm of spatial velocity vector v = y/U 2 + V 2 + W 2 for stars 
in our sample is shown on Fig. [6j The average value of this norm is 214 km/s with the 
standard deviation 86 km/s. Fig. [7] shows the V versus [m/H] dependence for single (left 
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upper panel) and double and multiple (right upper panel) stars in our sample. Middle and 
lower panels of this figure show the dependence of the multiplicity of our stars on four 
kinematic parameters UVW and v. The figure shows that the frequency of double and 
multiple stars in our sample weakly depends on the change of UVW components. We can 
draw lines corresponding to a constant multiplicity of approximately 35% within a 70% 
confidence level (see middle panel and left lower panel). For norm of spatial velocity vector 
a trend of decreasing multiplicity with increasing v is seen (see right lower panel). The 
data in the 95% confidence level do not contradict the hypothesis f sys tems(v) = const or a 
small growth of multiplicity with the increase of v. Nevertheless, with the probability of 
at least 70% we can say that with the increase of the norm of spatial velocity vector the 
frequency of double and multiple systems in our sample falls. Most likely, the bigger the 
spatial velocity of a star, the less the probability that it has a companion. 

Forty five from 223 stars in our sample are moving on retrograde galactic orbits 
(V < —220 km/s), i.e. in the opposite direction to the rotation of our Galaxy. For these 
stars, the ratio of single : binary: triple:quadruple systems is 32:12:1:0 and their multiplicity is 



Carney et al. 



(l2005aj ) . analyzing a sample of 374 stars on highly retrograde orbits 
(V < —300 km/s) showed, that the frequency of spectroscopic binaries among them is two 
times smaller than that for the stars moving along with the Galaxy's rotation. Our data, 
bearing an analysis of a wide range of periods, do not contradict their findings, yet, an 
insignificant number of highly retrograde objects in our sample does not allow us to make 
any definitive conclusions. 



7. Period distribution 



The distribution of orbital periods for binary stars in our sample is sho wn in Fig. [8] 
The periods of spectroscopic pairs are taken from iGoldberg et all (120021 ) and lLatham et al 
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( 120021 ). The periods of astrometric pairs were derived with the help of the generalized 
Kepler's third law on the basis of an empirical relation of the projected angular separation 
between the components and the semi-major axis. Knowing the system's parallax tt and 
the projected angular separation between the comp onents p, th e expe cted value of the 



semi-major axis is calculated using the formula from 



Allen et al 



f l2000l ) 



(a) = antilog 



log- + 0.146 

7T 



(4) 



where (a) is expressed in astronomical units, p and tc in arcseconds. Taking each system 
individually, we derived the sum mass of the components from the temperatures of 
the primary (CLLA) and secondary components, and from the magnitude difference, 
if the temperature of the secondary was unknown. To do this, we used models from 



Baraffe et al. 



( 119971 ). Angular distances p between the components of a strometric p airs 



wer e obtained from speckle inte r 



and 



Zapatero Osorio fc Martin 



erom etric observations or adopted from 



Allen et al. 



(120001 ) 



( 120041 ) . If the systems parallaxes were known from the 
HIPPARCOS catalog with an accuracy of better than 30%, we used them instead of the 
distances cited in the CLLA catalog. As a result, we were able to determine the periods for 
60 binary systems out of 64 in our sample. The periods for the four remaining suspected 
binary systems, — G186-26 and G210-33 and two bl ue stragglers, BP +25° 1981 and G43-3 



( iCarney et al. 



20011 ). are too long to be determined (lLathamll2008l ). The period distribution 
for 60 binaries and 10 multiple systems (18 subsystems of 9 triple stars and 3 subsystems of 
a quadruple star G89-14), is shown in Fig. [9j The distributions corrected for Opik effect 
and unresolved components on Fig. [8] and |9] are marked by a solid line. 

Let's compare the resulted distribution with an analogous one for the thin disc stars 
(Fig. ITOl) . The maximum of an unsymmetrical period distribution for the stars in our 
sample lays in the range of logP = 2 — 3 dex (i.e. hundreds of days). For Population I stars, 
the logP distribution, which Duquennoy and Mayor approximated by a Gaussian, has a 
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maximum in the range of 4 — 5 dex (tens of thousands of days). An important feature of our 
distribution is a small number of short period (logP < 1) p airs. This range is rep resented 



20021 ). which is 



by the only SB2 star G183-9 with a period of about 6 days ( iGoldberg et al 
excluded while accounting for the Opik effect. 

b make a ho mogeneity check of the two samples, we conducted a nonparametric x 2 



test (IKremer 



20071 ). By an accidental coincidence, the number of compared periods equals 
to 81 for each population. The test shows that the hypothesis of uniformity (i.e. common 
general population of two samples) can be rejected with a more than 95% probability. 

Both distributions (Fig. [10]) may be distorted by the selection effects. However, the 
differences in both the shape of the distributions and the location of the maxima could be 
induced by the dynamical evolution, undergone by the Population II stars. For example, 
a small quantity of old systems with periods of more than ten thousand days may be 
the result of a dynamical evolution at the stage of Galaxy formation. According to a 
recent concept, a bigger part of the stellar halo in our Galaxy was formed from small 



galactic systems (IBell et al. 



20081 ). It is quite likely that at the stage of accretion of small 
galaxies-satellites of the Milky Way and their destruction, the physical conditions were 
favourable for dissociations of wide pairs with low binding energy. Attempting to explain 
the period distribution of Population II stars by a destructive impact of giant molecular 
clouds and other local perturbations of th e gravitational potential of our Galaxy on the old 



stars orbiting around the galactic centre (I Weinberg. Shapiro fc Wasserman 



19871 ). is quite 



problematic. Such objects spend most of their lifespans away from the galactic plane. 
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8. Interesting triple and quadruple systems in our sample 



In this section we will examine two old multiple systems that we find remarkable: 
a triple system G40-14 and a quadruple G89-14. Such objects are of great interest for 
the studies of the dynamical evolution and for the checks of various criteria of dynamical 
stability. In Table HJ we listed all the detected systems from our sample having more than 
2 components. From ten multiple systems, nine are triples and only one is quadruple. 

The uniqueness of the triple system G40-14 is in its retrograde galactic orbit 
with V ~ —230 km/s (CLLA). The inner subsystem of G40-14 is an SB1 pair with a 



period of 60.615 days ( ILatham et al. 



2002). The outer subsystem is formed by a visual 



component, which is located 98" away from the spectroscopic pair. Assuming that the 



system's he 



( jAllen et al. 



iocentric distance is 235 pc, the expected semi-major axis is (a) ~ 30 000 AU 



20001 ) . We made a check for re trograde o 



August 2007) of the Multiple Star Catalog ( iTokovinin 



rjects in the latest version (dated 13 



19971 ). which is a compilation of 1158 



known stellar systems with three or more components. In order to do that, we calculated 



the U V W components of spatial velocities for the catalog s tars from the cited par a 



radial velocities and proper motions, using the formulae from I Johnson & Soderbloml ( 119871 ). 



laxes, 



It appeared that only one object, a triple system ADS 16644, is moving on a retrograde 
galactic orbit (V ~ —330 km/s). Hence, G40-14 is the second of all known systems with 
more than two components moving against the rotation of the Galaxy. 



G89-14 (for details, see 



Rastegaev 



20091 ) is a system with the highest multiplicity in 



our s ample. It consists of four components (Fig. [TTj) : an SB1 pair AB with a period of 190 



days (ILatham et al 



this SB1 pai r ( 



(lAllen et al. 



200 2}), a speckle interferometric component C located at ~ 1" from 



Rastegaev et al. 



20071 ) and a common proper m otion companion D at 34" 



2000). Based on the data from 



Allen et al. 



from iBaraffe et al 



( 120001 ) . the evolutionary tracks 



( 119971 ). and on our speckle interferometric measurements, we evaluated 
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the period ratio of the three G89-14 subsystems: 0. 52 : 3 000 : 65 000 yr. Another 
well-known metal-poor quadruple system is NQ Ser (jTokovininl 



(INordstrom et al 



19971 ) with [Fe/H] = -1.05 



20041 ). This multiple star was r epeatedly observed on the BTA by means 



Balega et al. 



20061 ) . According to our information, G89-14 



of speckle interferometry (e.g., 
with [m/H] = —1.9 (CLLA) is the most metal-poor quadruple system known to date, which 



makes it an interesting object for a more detailed study. 



9. Multiplicity of metal-poor stellar streams 



Stellar streams (e.g., lEggenl ll996al Jbl) are associations of stars possessing similar 
kinematics and metallicity. The study of such streams allows restoring to a certain degree 
the picture of the formation of various dynamical structures in our Galaxy. Traditionally, 
the stellar streams are being selected in a certain phase space and then their origin is 
interpreted using the data of spectroscopic analysis. In a phase space, a fine structure like 
stellar multiplicity can give additional information on the dynamical evolution of the stream 
and its primogenitor. However, until now it was not taken into due consideration. 



The following six stars of our sample: G1 0-4, G 



are part of the Kapteyn's star moving group (jEggen 



3-9, G 60-48, G24-3, G18-54, G28-43, 



1996ah . 10 other objects: G130-65, 



G75-56, G5-35, G40-14, G11 4-25, Gll-44 , G13-35, G183-11, G182-32, G126-52, belong to 



the Ross 451 moving group (jEggen 



1996bl ). In Table Owe are listing some characteristics 
of these two halo streams. In the penultimate column of the table, you can see the ratio of 
single, binary and triple systems for the group members in our sample. In the last column, 
an analogous estimate is given for all known members of the groups, the data taken from 
literature. Unfortunately, multiplicity of these stars is poorly studied and the ratio in the 
last column can only serve as a lower limit for the frequency of binary and multiple systems 
in the streams. The two moving groups listed above have comparable multiplicities, both 
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exceeding 10%. 



The fact that we find similar multiplicities in moving groups, does not contradict the 
dynamical hypothesis of their origin. In this case, the stellar streams are formed by a 
random selection from the general population of field sta rs. Recent works on the problem of 



origin of stellar streams show that the Herc ules stream (IBensbv et al. 



Pleia des, Hyades, and Sirius moving groups ( IFamaey et al 



2007 



2007), as well as the 



Famaey. Siebert &: Jorissen 



20081 ). formed as a result of dynamical (resonant) influence of our Galaxy on the field stars. 
However, the scenario of accreted stellar streams cannot be ruled out. This requires similar 
multiplicities of the progenitors of the flows. Further detailed studies are required to help 
answer the question whether there exist any distinctions between various stellar streams. 



10. Conclusion 



In this paper we examine a sample of 223 subdwarfs belonging to the F, G and early K 
spectral classes, located within 250 pc from the Sun, with the metallicities [m/H] < — 1 and 
proper motions // > 0.2"/year. Stars make up about 5% of the total number of objects of 
this type in the studied space volume. The subdwarfs were o bserved using the spec t rosco pic 



([Goldberg et al 



2002 



Latham et al 



2002 ), interferq metric ( IRastegaev et al. 



and visual methods (jZapatero Osorio &: Martin 



2007 



20081 ) 



20041 ). Presented sample is most thoroughly 



studied in terms of stellar multiplicity in a wide range of orbital periods (orbital axes) 
among the Population II field stars. 

As a result of observations of the sample stars using the method of high angular 
resolution on the BTA, we detected 20 speckle interferometric components for 19 primaries. 
Seven of them were known as spectroscopic pairs and four — as astrometric binaries (Fig. 
HJ). Nine systems were resolved for the first time: 5 binaries (G191-55, G114-25, G142-44, 
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G28-43, G130-7), 3 triples (G87-47, Glll-38, G190-10) and one quadruple (G89-14). 



Combining different research methods allows us to estimate the frequency of the 
systems of different multiplicity with the orbital semi-major axes ranging from a few to 
tens of thousands of astronomical units. The ratio of single, binary, triple and quadruple 
systems among 221 primary components in our sample amounts to 147:64:9:1. More than 
half of the stars in the sample are members of binary and multiple systems: 

2B + 3T + 4Q 



f stars 

Multiplicity of the sample is 



S + 2B + 3T + 4Q 



52% ± 6%. 



B + T + Q 



systems 



33% +6% 



7%- 



S+B+T+Q 

As before, a 95% confidence interval was used as the error of obtained multiplicity. For 
spectroscopic systems, we made an analysis of the number of undetected components using 
analytical calculations, while the similar estimates for astrometric pairs were obtained 
using the Monte Carlo numerical method. For our sample the corrected for undetected 
components and selection effe cts (eventually only the Op ik effect was accounted for) 



multiplicity is f syst e 
disc G dwarfs. 



40%. 



Duquennoy fc Mayor! (Il99ll ) give f systems ~ 60% for the thin 



For colder subdwarfs of K-M spectral classes, IJao et al.l (120091 ) deduced 
S:B:T:Q=46:12:2:2, and f systems = 26 ± 6% accounting for spectroscopic, speckle in- 



ter 
by 



erometric and visual data. Within errors our result coincides with the result obtained 



Jao et al. 



(120091 ) . 



Seven stars in our sample are blue stragglers. We found that with an exception of 
G245-32, all of them are binaries. This supports the hypothesis of the connection of the 
blue stragglers phenomenon with their binary nature. 
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The ratio of binary:triple:quadruple systems (B:T:Q) among the Population II stars 
in our s ample is 64:9:1. This can 



I stars ( IDuquennoy fc Mayor 



De compared with the ratio 40:7:2 for the Population 
19911 ) . The difference between the ratios is statistically 
indistinguishable. Therefore, we came to a conclusion that a stable hierarchical multiple 
system is the universal evolutional outcome of the star formation process both at the time 
of the creation of our Galaxy and nowadays. 

Three of the resolved by us binary systems, G76-21 (HIP 12529), G114-25 (HIP 44111) 
and G217-8 (HIP 115704), have very low metallicities ([m/H] < —2). A triple system 
G40-14 also belongs to this metallicity range (see Table H]). Altogether, 63 primaries from 
our sample have [m/H] < —2. From these, 18 primaries have one companion and one, 
G40-14, has two of them. To date, only a few high multiplicity (N > 2) systems in very low 
metallicity regime are known. Further accumulation of empirical data for these objects will 
help answer the question about a possible dependence of the orbital elements distribution 
of binary and multiple systems on their metallicity. 

We did not find any significant differences in the multiplicity r atios of sub dwarfs 



moving on prograde and retrograde (V < —220 km/s) galactic orbits. ICarney et al.l (J2005aJ) 
found a decreased ratio of strongly retrograde (V < —300 km/s) binaries: 10% ±2% against 
28% ± 3% for a prograde sample. Their conclusion is supported by our study of 11 stars 
with V < —300 km/s: only one of them is a binary. With the increase of the norm of 
spatial velocity vector v = \/U 2 + V 2 + W 2 the frequency of double and multiple systems 
in our sample falls with a probability of at least 70%. 

We have shown that the distribution of orbital periods of the old stars differs both in 
shape and in the maximum's location from that of Population I stars (Fig. ITOl) . Most of the 
detected Population II binaries have periods between 1 and 10 yrs. The period distribution 
for G, K and M thin disc dwarfs does not depend on a spectral class. In the semilog scale, 
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it can be approx imated by a Gaussian with the maximum at logP„ 



in 



Kroupa 



5 dex (see Fig. 1 



19951 ). Orbital periods of the thin disc stars are distributed more symmetrically. 
Compared to Population II stars, the maximum of the distribution is shifted towards larger 
P by two orders. 

An important feature in the period distribution of binary and multiple subdwarfs is 
the lack of short-period systems with the periods of less than 1 day. The range of P < 10 
days is represented in our sample just by one SB2 system — G183-9. The reason for that 
lack of old short-period pairs could be in the dynamical evolution of the systems with short 
orbital periods, whic h could lead to a merge of the components and to the formation of 



blue stragglers (e.g., 



Bailyn 



19951 ) . Another important peculiarity in period distribution 
of old metal-poor stars is the presence of couples with P > 10 s days (Fig. [8] and [9]). It is 
quite possible that such enormous periods could be an error occuring from the formula (jlj). 
Such low-binding energy objects survived over many billions of years and bear important 
information on the mass density distribution in our Galaxy. The orbital periods of the 
subdwarfs and thin disc dwarfs lie nearly in the same range and constitute 10 orders. It is 
impossible to explain such a wide range of periods by the dynamical evolution only, as it is 
forme d at the earliest stages o 



(e.g., 



Kroupa fc Burkert 



the stellar system formation in the nuclei of molecular clouds 
20011 ) . Our data show that the range of possible orbital periods 
of binary and multiple systems is not decreasing over billions of years of the dynamical 
evolution. Only the shape of the period distribution is changing. An interpretation of the 
distinctions between the period distributions of the stars of different populations requires 
further study. 

The substellar mass companions (bro wn dwarfs and planets) s tayed beyond the scope 



of our consideration. Current studies (e.g., 



Fischer fc Valenti 



20051 ) testify to the existence 



of a correlation between the metallicity of stars and the presence of orbiting planets. At the 



-27- 



time of writing, catalogs of stars with planets (lhttp://exoplanet.eul ; lhttp://exoplanets.org| ) 
do not contain stars with [Fe/H] < —1. Some researchers claim that brown dwarfs and 
some very low mass stars (and possibly planets?) , form a separate p opulation with its own 



Kroupa et al. 



20031 ). As likely as not, the 



multiplicity and kinematical properties (see, e.g. 
low metallicity regime may influence the formation of this population. It is quite possible 
that Population II field stars do not contain any brown dwarfs or planets as companions at 
all. Another important question is whether there exist any substellar components in the 
halo and thick disc systems of high multiplicity (N > 2). 

The ratio of binary and multiple systems among the Population II stars, found in this 
study, does not contradict the hypothesis that the chemical composition of protostellar 
molecular clouds makes but an insignificant impact on the star formation process. This 
indicates that the halo stars were formed as a result of fragmentation of molecular clouds' 
nuclei, similarly to the way the stars form today. However, this issue remains unclear when 
we consider the substellar mass regime. 

The frequency of binary and multiple metal-poor stars imposes some restrictions on 
the formation of the stellar halo in our Galaxy as well. There are two general scenarios of 



the formation of the Milky Way's stellar halo (see lMajewskilll993l for more details): 



Most of the halo stars were born in globular clusters or dwarf galaxies, which 
were then accret ed and destroyed in the gravitational potential of the Milky Way 



f lBell et al 



20081 ) 



The halo stars are genetically bound with our Galaxy. The accretion of globular 



clusters and dwarf galaxies observed today ( jlbata. Gilmore fc Irwin 



19951 ). produces 



only a small fraction of halo stars. The 



stellar halo ( IHelmi et al. 



1999 



dnematical structures detected in the 



Bell et al 



20081 ) are a consequence of the Galaxy's 



gravitational potential inhomogeneities and a manifestation of various resonances. 
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Any scenario of the stellar halo formation has to impose certain restrictions on the rate 
of binary and multiple systems and on their characteristics, i.e. on the distributions of 
orbital periods, component mass ratios, eccentricities, etc. Particularly, high percentage 
of binary and multiple halo field stars indicates that the formation of stellar halo via 
the destruction of globular clus ters is unlikely in ou r Galaxy, since the relative number 



of binaries in globular clusters (ISollima et al. 



20071 ) is smaller than that of metal-poor 



field stars. It is currently not c 



influe nce the binary frequency (jlvanova et al 



ear how does 



;he dynamical evolution of globular clusters 



2005 



Hurley. Aarseth & Shara 



2007 



Sollima 



20081 ). therefore the scenario of the halo field subdwarfs formation through a dissociation of 



globulars cannot be fully discarded. 

The question of the differences between stellar streams in terms of binary and multiple 
systems requires further accumulation of observational data. Our material does not 



contradict neither an assumption of the parity of binary and multip 
different streams, nor the hypothesis of their dynamical origins (e.g., 



e stars frequency in 



Famaey et al 



20081 ). 



Some of the detected speckle interferometric pairs, G76-21, G63-46, G28-43, G217-8, 
G130-7, G102-20, BD+19° 1185A, G87-47, with presumably short orbital periods are 
suitable for monitoring for orbit calculations and mass determination of the metal-poor 
stars. These studies can contribute to a calibration of the mass-luminosity relation and to 
the verification of the theories of dynamical evolution. 
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A. Impact of Galactic structure on spatial distribution of sample stars in 

solar neighbourhood 

To be able to estimate the number of stars in a given volume by an extrapolation of the 
calculated quantity of stars in a smaller volume, we have to take into account the structure 
of the Galaxy. A null hypothesis is that in the volume of 250 pc from the Sun, the amount 
of thick disc stars is 1 000 times bigger than that within 25 pc. We assume a uniform 
distribution of halo stars on the scale of 100 pc 3 in the solar vicinity. Let us consider how 
an exponential decrease in the number density with an increase in the distance from the 
galactic plane in the sample of thick disc stars affects the null hypothesis. 

Let's introduce a k = Ni/N Q coefficient, where 

270 

Ni = J (250 2 - (z - 20) 2 ) x e-^ /H dz. (Al) 

-230 

The Ni is proportional to the quantity of the thick disc stars within 250 pc from the 
Sun. We assume that the S un is located at the distance of 20 pc above the galactic plane 



( [Humphreys fc Larsen 



1995T ). The quantity of stars, located within 25 pc from the Sun is 



proportional to: 

45 



N = J (25 2 - (z - 20) 2 ) x e- ]z]/H dz. (A2) 



-5 



A true quantity of the stars in these volumes can be obtained by a multiplication of A^i 
or N on tt x p , where po is the number density of thick disc stars at the galactic plane 
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(z — 0). Apparently, in case of uniform distribution of stars in space, k is equal to 1 000. 
All values in these formulae are expressed in parsecs. The location above the Galactic plane 
is designated by z. The thick disc scale height H va ries from 1000 to 1500 pc. For our 



calculations we took H = 1048 pc (jVeltz et al. 



20081 ) . After integrating the right parts of 



(I All) and ( 1A2I) . we obtain k = 933. Therefore, the deviation from the homogeneous number 
density distribution of thick disc stars within 250 pc from the Sun is approximately 7%. 
The median metallicity of the stars in our sample is [m/H] w —1.6. The percentage of the 
metal-weak thick disc tail stars in t he range of —1.6 < [Fe/H] < — 1 may reach 60 — 70% 



( [Morrison. Flynn fc Freeman 



1990 



Beers Sz Sommer-Larsen 



19951 ). Nevertheless, while 



considering our sample, the Galactic structure in the solar neighbourhood can be neglected. 
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Fig. 1. — Metallicity [m/H] versus the ^-component of spatial velocity for the stars in our 
sample. The circles represent the halo stars and the squares mark the thick disc stars. The 
dashed line separates the stars on prograde (upper half) and retrograde (lower half) orbits. 
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Fig. 2. — (S— V) - My diagram for 213 stars in our sample (crosses) . To make a comp arison 



we plotted the stars from the M13 globular cluster (d ots) from 



Rev et al. 



distance modulus for M13 is (m — M) = 14.38 ± 0.10 (IGrundahl et al 
stragglers from our sample are marked with filled circles. 



(120011 ). The 



19981 ). Seven blue 
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Fig. 3. — A compar ison between the trigonometric distances from the HIPPARCOS catalog 



(Ivan Leeuwen 



20071 ) and photometric distances from the CLLA for 133 sample stars. The 
dotted line is a linear regression y = —18.58 + 2.23a;. r = 0.42 is the correlation coefficient. 



The solid line represents the y = x function. 
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Fig. 4. — Schematic representation of the companions in our sample detected using different 
methods. 
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Fig. 5. — Multiplicity versus metallicity for stars in our sample (solid line). Dashed lines 
mark the 95 % confidence interval. Dotted lines mark the 70 % confidence interval. 
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Fig. 6. — Norm of velocity vector distribution for stars in our sample. 
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Fig. 7. — The upper panel shows V versus [m/H] dependence for single (left part) and double 
and multiple (right part) stars. Middle and lower panels show the frequency of double and 
multiple stars depending on the components of the velocity vector and its norm. Dashed 
lines mark the 95% confidence interval. Dotted lines mark the 70% confidence interval. 
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Fig. 8. — Period distribution for 60 binary stars in our sample (dashed line). The solid line 
marks the same distribution but corrected for Opik effect and unresolved components. Error 
bars represent square roots from the number of stars in each bin. 
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Fig. 9. — Period distribution for 60 binary stars and 21 subsystems of 10 multiple stars in 
our sample (dashed line). The solid line marks the same distribution but corrected for Opik 
effect and unresolved components. Error bars represent square roots from the number of 
stars in each bin. 
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Fig. 10. — Distribution of densities of orbital periods for binary stars in our sample (Pop- 
ulation II) interpolated by a c ubic spline and idem f or the thin disc stars (Population I), 
approximated by a Gaussian in iDuquennoy fc Mayor! (119911 ) . 
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Fig. 11.— On the right: CPM subsystem of G89-14 (POSS archive). On the left: the speckle 
interferometric subsystem (image reconstructed from speckle interferograms using bispectral 
technique). The magnitude difference between AB and C is 4.2 m in 800/100 filter. 
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Table 1: Properties of three stellar Populations 



Value 


Thin Disc 


Thick Disc 


Halo 


&u 


43 


67 


131 


(V) 


-15 


-53 


-226 


°v 


28 


51 


106 


a w 


17 


42 


85 


([Fe/H]) 


-0.1 


-0.8 


-1.8 


& [Fe/H] 


0.2 


0.3 


0.5 


Fraction/ 


0.925 


0.070 


0.005 



According to the data from Robin et al. (2003). 
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Table 2: Multiplicity function of metallicity 



[m/H] 


(-3,-2.5] 


(-2.5,-2.0] 


(-2.0,-1.5] 


(-1.5,-1.0) 


S : B : T : Q 

J systems 

Number of systems 


17:5:1:0 

26% +23% 

ZO/0 -16% 

23 


26:13:0:0 

oo % +17% 
<J^/o_ 14% 

39 


41:20:2:1 

36% +13% 

— 12% 

64 


62:26:6:0 
94 



* 95% confidence interval is given 



Table 3. Multiplicity as a function of kinematics 



U (km/s) 


(-300, -200] 


(-200,-100] 


(-100,0] 


(0, 100] 


(100, 200] 


(200, 300] 


S : B : T : Q 

j systems 

Number of systems 


11:5:1:0 
35%+ 27% 
17 


29:14:2:1 

3 ' /c -14% 
46 


40:14:1:0 

27/c -ll% 

55 


32:17:5:0 

41/o -13% 
54 


24:8:0:0 

ncCK + 18% 
^0/0_ 14% 

32 


9:4:0:0 

31%+-% 

13 


V (km/s) 


(-300, -250] 


(-250,-200] 


(-200, -150] 


(-150, -100] 


(-100, -50] 


(-50,0] 


S : B : T : Q 

f systems 

Number of systems 


9:6:0:0 

40/0 -24% 
15 


25:8:1:0 

9fi%+ 18% 
26/c_ 13% 

34 


34:15:1:0 

32/0 -12% 

50 


27:10:4:0 

34/0 -14% 
41 


25:11:1:1 

34/0 -14% 

38 


13:10:1:0 

4fi qz,+21% 
46/0 -20% 

24 


W (km/s) 


(—300, —100] 


(—100, —50] 


(—50, 0] 


(0, 50] 


(50, 100] 




S: B:T :Q 

/systems 

Number of systems 


11:2:1:0 
21/o_ 16% 
14 


18:13:1:0 

44% +18% 

32 


41:18:4:0 

35% +13% 
63 


40:13:2:1 

29% +13% 

56 


29:14:1:0 

•34/° -14% 
44 




VU 2 + V 2 + W 2 (km/s) 


(0, 100] 


(100, 200] 


(200, 300] 


(300, 400] 


> 400 




S: B:T :Q 

f systems 

Number of systems 


10:8:2:0 
50%^3% 

20 


48:23:5:1 
77 


57:24:2:0 

31%+ ll0/f ' 
,31/0 -9% 

83 


26:7:0:0 
21% +18% 
33 


5:0:0:0 

Q% +52% 
U/0 -0% 

5 





Note. — For /systems the 95% confidence interval is given. 
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Table 4: Population II multiple systems in our survey 



Name 




[m/H]* 


Distance* 

(PC) 


Multiplicity 


G95-57 


8 m 78 


-1.05 


25 


3 


BD+19° 1185 


9.3 


-1.47 


42 


3 


G89-14 


10.4 


-1.9 


94 


4 


G87-45 


11.44 


-1.49 


123 


3 


G87-47 


10.34 


-1.34 


62 


3 


Glll-38 


9.11 


-1.04 


50 


3 


G40-14 


11.2 


-2.71 


164 


3 


G59-1 


9.52 


-1.14 


50 


3 


G17-25 


9.63 


-1.54 


35 


3 


G190-10 


11.22 


-1.92 


91 


3 



* According to the data from CLLA. 
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Table 5: Moving groups in our sample 



Group 


[Fe/H] 


U 


V 


W 


n x : n 2 : n 3 


Ni : N 2 : ^ 3 


Ross 451 


-1.5 


-89 ± 104 


-346 ± 5 


-37 ± 77 


8:1:1 


36:5:1 


Kapteyn's 


-1.5 


13 ±82 


-288.5 ±6.0 


-16 ±67 


4:2:0 


29:4:0 


star 















Data are taken from Eggen (1996a,b). 



